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Abstract 

In this paper, we study long-term correlations and multifractal properties elaborated from time 
series of three-phase current signals coming from an industrial electric arc furnace plant. Implicit 
sinusoidal trends are suitably detected by considering the scaling of the fluctuation functions. Time 
series are then filtered via a Fourier-based analysis, removing hence such strong periodicities. In 
the filtered time series we detected long-term, positive correlations. The presence of positive cor¬ 
relations is in agreement with the typical V-I characteristic (hysteresis) of the electric arc furnace, 
providing thus a sound physical justification for the memory effects found in the current time series. 
The multifractal signature is strong enough in the filtered time series to be effectively classified as 
multifractal. 

Keywords — Multifractal analysis; Electric arc furnace; Time series analysis; Fourier detrending. 


1 Introduction 

The analysis of complex systems 13 [HI HO] is recently gaining traction with the development of new 
computational and theoretical tools, specifically designed to tackle the problem of disentangling the 
complex interactions among different parts of a process. A common approach for this purpose is the 
study of scale-invariant properties |37| . i.e., those features that manifest in a similar manner at all 
magnification levels of the system. The scale-invariance property is a peculiar feature of fractal structures, 
and it is quantifiable by evaluating a suitable collection of scaling exponents that compactly describe 
how the system evolves when observed at different scales m- In particular, the scaling of the correlation 
properties of a stochastic process could provide valuable information on the underlying dynamical system, 
especially in cases in which the observed system denotes long-term memory and nonlinear behaviors 

[31IH]- 

Long-term memory properties of a stochastic process are firstly described by power-law like decreasing 
of the autocorrelation calculated for suitable observables. This feature can be characterized by the 
value of the so-called Hurst exponent H [35], which assumes values in [0,1]. The value of H is 0.5 
when the process corresponds to uncorrelated noise, whereas if the process is persistent (correlated) or 
antipersistent (anticorrelated) then H is greater than or less than 0.5, respectively. However, conventional 
methods employed to analyze the long-term correlations of a time series (e.g., spectral analysis and 
Hurst analysis) are known to be misleading when data are non-stationary and contains deterministic 
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trends 13 [E]. Therefore, in many experiments it is important to distinguish intrinsic fluctuations 
characterizing the process from pure trending behaviors. A well-established approach employed to this 
end is the Detrended Fluctuation Analysis (DFA) |36], which has been successively generalized in the 
so-called Multifractal Detrended Fluctuation Analysis (MF-DFA) [EllinilSH]- MF-DFA accounts for the 
possibility to characterize different degrees of correlation by considering a collection of scaling exponents. 
Such scaling exponents provide a way to perform classical multifractal analyses on time series. Notably, 
it is possible to calculate the so-called multifractal spetrum, which describes all relevant multifractal 
characteristics of the system under analysis. MF-DFA technique has been used in a multitude of contexts, 
such as breast cancer imaging [12j , river flow modeling [26] , brain signals [10] , solar filaments and sunspots 
[HISS], financial returns [13 m], high-energy physics [13 , earth’s gravity [33 j protein contact networks 
[13, and protein molecular dynamics [33- 

Detrending is a very crucial aspect in multifractal analysis of time series [IIEIIIT]. In fact, trends 
are usually a source of spurious multifractal signatures observed in the analyzed series. Conventional 
MF-DFA is able to deal with trends compatible with (higher-order) polynomials. Nonetheless, trending 
could appear also as a consequence of periodicities, such as those induced by seasonality for long time 
series. To this end, Chianca et al. [3 proposed to perform a preprocessing stage involving a Fourier 
analysis in order to identify and then remove specific frequencies with large amplitude. Such a filtering 
technique demonstrated to provide good results in several different contexts [23 dZl 113 > while however 
criticisms and limitations of the method have been highlighted in the context of sunspot time series [14] . 

In this paper, we study long-term correlations and multifractal properties elaborated from time series 
of three-phase current signals coming from an electric arc furnace (EAF) [23111] • EAF modeling and 
analysis of related data is an active area of research [3[I71I23IM3- Previous studies focused mostly on 
system identihcation and prediction problems with the aim of power optimization and control. However, 
to our knowledge there is no track in the literature of studies employing MF-DEA to characterize EAF 
signals. An EAF is a plant utilized in industrial environments for the production of steel by a melting 
process of metallic scraps. EAF consists of a cylindrical melting pot covered by a removable vault, inside 
of which are placed one or more electrodes. The scraps to be molten (metallic charge) are deposited on 
the bottom of the melting pot and, by applying AC or DC tension to the electrodes, an electric arc is 
generated between the extremities of such electrodes and the scraps. The current traversing the metallic 
charge raises its temperature until the fusion point is reached, initiating the melting process. The V-I 
characteristic of the EAF is the result of a nonlinear process with memory m, since the electrical current 
value is determined as a function of both the voltage and the internal state, which in turn is determined 
by the past states of the system. Such a process gives rise to an electrical phenomenon called hysteresis 
0 . 

Toward removing all implicit periodicities proper of an alternating current setting, the EAF time 
series are preprocessed via a Fourier-based analysis in order to remove the first few complex coefficients. 
By analyzing the filtered time series we detected long-term, positive correlations. This outcome is in 
agreement with the typical V-I characteristic of the EAF, which shows a clear hysteretic behavior, as 
in fact the equivalent dipole of the arc is a dynamic and nonlinear one [3. The multifractal signature 
is consistent among the three-phases and it is strong enough in the filtered time series to be effectively 
classified as multifractal. Finally, we show that shuffling the series reveals that the multifractal signature 
is entirely due to the long-term correlations. 

The paper is structured as follows. Section [3 provides an introduction to the EAF context. In Section 
[3 we discuss our results regarding long-term correlations and multifractal properties of the analyzed 
three-phase current signals. Section [^provides the concluding remarks and future directions. Finally, in 
Appendix [^ we introduce the MF-DFA technique. 


2 Electric arc furnace and the V—I characteristic 

The analyzed EAF [31] is an alternating current furnace powered by a three-phase electrical supply. The 
plant contains three electrodes placed on the three vertices of a triangle, supported by vertical guides 
that allow controlling their distance with respect to the underlying metallic charge. When the applied 
voltage is sufficiently high to break the air resistance, an electric arc is struck from the extremity of each 
electrode to the metallic charge. In an ideal setting, the electric arc is created each time the absolute 
value of the alternating voltage exceeds a threshold and, consequently, the current assumes an almost 
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sinusoidal waveform. However, in real-world situations, the V-I characteristic is subject to many sources 
of noise, such as arc instabilities, non-monochromaticity of the input signal, and metallic charge cave-ins 
that hinder the analysis of the measured data. The latter effect, for example, might cause sudden drops 
in the current magnitude. 

The V-I characteristic of the system forms a closed trajectory in the tension-current phase space (see 
Fig. 0 . This behavior denotes a branching process in the joint evolution of those two quantities and 
it is the distinctive feature of a system subject to hysteresis [S]. This implies that the system possesses 
memory of the previous states. As we will show in this study, the hysteresis can be investigated through 
the analysis of the long-term correlations of the related EAF current time series. 



Figure 1: Hysteresis of the V-I characteristic of EAF. Image courtesy of Sadeghian and Lavers 


3 Results and discussion 

Our data consists in three current time series, each related to one of the three phases of the EAF under 
analysis m- In the following, we will denote the three phases as A, B, and C. Each time series is 
composed 1150000 samples, having a sampling frequency of 1920 sample per second. The inter-sample 
time is hence 520.833 x 10“® seconds, which gives a total of roughly 10 minutes of recording. All data 
are standardized prior to processing. 

Fig. i shows the autocorrelation |2(a)| and power spectrum |2(b)| for phase A. Autocorrelation is in 
agreement with a typical sinusoidal waveform of a current signal (60 Hz). The power spectrum offers an 
insight on the fact that such a signal contains many harmonics of significant power, which correspond to 
strong periodicities in the series. This is fairly expected, due to the nature of the physical system under 
analysis. However, in order to avoid the detection of spurious multifractal signatures, it is necessary 
to eliminate such periodicities. The interest in the analysis of such EAF signals lies hence in the non- 
linearities that are super-imposed during the melting process on the baseline sinusoidal signal. 

In Fig. [^we show the scaling of the fluctuations ([^ for the three original signals. A cross-over is visible 
in Fig. I^at a scale of roughly 32-40 time instants, which corresponds (considering the aforementioned 
inter-sample time) to a regime of 48-60 Hz. The cross-over is consistent among the three phases, as 
in fact in the figure the three fluctuation functions overlap almost perfectly. Therefore, fluctuations, 
when calculated on the original signals, immediately pick up the main sinusoid of the signal due to 
the frequency of the electrical current. In order to remove traces of all periods in the three current 
signals, we applied Fourier detrending 0 EH EH ED]. The procedure consists in transforming the signal 
in the Fourier space with the fast Fourier transform. Complex-valued Fourier coefficients are then order 
according to their amplitude (from larger to smaller). In order to obtain a suitable signal to be processed 
with MF-DFA, we have zeroed the first 3000 coefficients (corresponding to roughly 1.5 seconds) and then 
antitransformed the signal back to the time domain. The result of the filtering can be appreciated in 
Fig. El In the figure, we show an excerpt of the original signals together with the identified trends (left 
panels), while in the right panels we report the corresponding filtered versions. It is worth noting that, 
as expected, the identified trend is not compatible with a single sinusoid, as in fact the trend in the signal 
is a combination of the main sinusoid with several harmonics of different power (see Fig. |2(b)[). As we 
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(a) Autocorrelation function. 


(b) Power spectrum. 


Figure 2: Autocorrelation [2(a)| and power spectrum |2(b)| of current time series of phase A. It is possible 
to recognize the presence of several harmonics with significant power. 
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Figure 3: (colors online) Scaling of fluctuations of original (standardized) time series. Fluctuations, 
calculated with a linear fitting and q = 2, denote a cross-over compatible with the 60 Hz of the electrical 
current. 


will demonstrate in the following, such filtered signals contain long-term, positive correlations and a fair 
multifractal signature. 

Filtered signals are then analyzed with the MF-DFA procedure, which is described in Appendix]^ 
MF-DFA has been executed by using scales from 16 up to 16384, q ranging in [-5, 5], and with a linear 
fitting for the local detrending; same outcome has been obtained with a quadratic fitting. Results are 
reported in Fig. Scaling of fluctuations, shown in Fig. |5(a)[ provides a confirmation that the filtering 
stage removed all relevant periodicities in the original signals. In fact, log-log plot of the hltered signals, 
for all three phases, is now largely consistent with a global power-law scaling. Let us now focus on 
Fig. |5(b)[ which shows the generalized Hurst exponents Q. All three signals denote long-term, positive 
correlations: H{2) ~ 0.75 for phase A, H{2) ~ 0.78 for phase B and C. In Fig. 5(c) we show the local 
slopes of the fluctuation functions for g = 2. Local slopes of the three signals fluctuate around the 
estimated Hurst exponents (the global slope of the entire fluctuation function). As typically found in 
synthetic time series |36j . higher deviations are observed when moving to larger scales. The discovery of 
long-term, positive correlations in all three (filtered) electrical current signals is well-justified and it is in 
agreement with the typical V-I characteristic of EAF (Sec. [^. In fact, hysteresis is a phenomenon that 
is observed when the output of a (nonlinear) system does not depend only on its current input, but also 
on past system states. Such a memory effect is imprinted in the EAF time series, as highlighted by the 
generalized Hurst exponents shown in Fig. |5(b)| 

Another important fact in multifractal analysis is the distinction between monofractal and multifractal 
scaling: if the scaling of a system can be fully characterized by a single exponent (Hurst exponent in 
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Figure 4: (colors online) First 500 samples of the time series (left panels, black lines indicate the original 
signal while in red we show the trend) and related filtered version (right panels, blue lines). Phase A, B, 
and C shown from top to bottom. 


our case) then it is monofractal, otherwise it is called multifractal. The analysis of the multifractal 
spectrum shown in Fig. |5(d)| pr ovides the answer to such a question. All three signals denote a fair 
multifractal scaling: the width (|lT]) of the spectrum is Aa ~ 0.30 for phase A, Aa ~ 0.21 for phase B, 
and Aa ~ 0.31 for phase C. It is worth mentioning for phase A and C a small but yet visible tendency of 
being insensitive to larger fluctuations, which can be deduced from the fact that the spectra are slightly 
right sided. The converse holds for phase B, which instead appear as slightly left sided - as in fact the 
corresponding Hurst exponents tends to be flat for negative q. As discussed by Drozdz and Oswi§cimka 
[8], right sided spectra are a less usual in experimental time series. The authors [8] offer a mechanism for 
reproducing such an effect by constructing a synthetic time series starting from a binomial cascade and 
successively super-imposing a white noise signal considering a suitable threshold for editing the original 
cascade. However, we believe that in our case such an asymmetry is likely to be a byproduct of the 
preliminary filtering stage performed via Fourier coefficient truncation. Our results indicate that all 
three (filtered) time series are compatible with signals having both long-term, positive correlations and 
a fair amount of multifractal features. This suggests that the non-linearities and memory effects of the 
melting process are super-imposed on the original current signals as noise. Such a noise retains relevant 
information regarding the physical process, which we suitably detected by means of MF-DFA. 

Multifractality could be observed as a byproduct of the long-range correlations and/or it could emerge 
as a consequence of an underlying broad probability distribution m- Shuffling the time series provides 
hence a reliable proof to disambiguate the origins of the observed multifractality. In fact, if the multifrac¬ 
tal behavior is entirely due to the long-term correlations, then the shuffled signals would not present any 
multifractality at all; shuffling the data destroys all long-term correlations. This can be quantitatively 
verified by observing a relevant shrinking of the multifractal spectrum width. On the other hand, if 
multifractality is (also) due to a broad underlying density, then shuffling would have only a visible but 
yet limited impact on the width of the resulting multifractal spectrum. As shown in the inset of Fig. 
|5(d)[ the shuffled signals do not contain any sign of multifractal signature. This leads us to claim that the 
observed multifractal properties of the current time series are entirely due to the long-term correlations 
- the memory effects due to the hysteresis. We argue that a broad probability density function, in our 
case, might not be possible, due to the EAF physical constraints. In fact, EAF operates in an industrial 
setting that is suitably controlled in order to achieve specific operating conditions. To summarize, in 
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Table we report all relevant coefficients determined via MF-DFA. 




(a) Scaling of fluctuations. 


(b) Generalized Hurst exponent. 



Scale 


(c) Local slopes of fluctuations for q = 2. 



(d) Multifractal spectrum. 


Figure 5: Scaling of fluctuations |5(a)[ generalized Hurst exponents |5(b)[ local slopes 5(c)| of the fluctuation 
functions for q = 2, and multifractal spectra |5(d)| of the filtered time series. In Fig. |5(a)[ fluctuation 
functions are shown for all considered q and vertically shifted to improve visualization. Horizontal lines 
in Fig. |5(c) I denote the respective estimated Hurst exponents of each phase. The inset in Fig. |5(d)| shows 
the spectra obtained after shuffling the time series. 


Table 1: Summary of all relevant exponents and coefficients found with MF-DFA. 



phase A 

phase B 

phase C 

H{2) 

0.75 

0.78 

0.78 

r(2) 

0.51 

0.56 

0.56 

Aa 

0.30 

0.21 

0.31 

Aa (shuffled) 

0.04 

0.06 

0.05 


4 Conclusions and future directions 

In this paper, we studied the long-term correlations and multifractal properties elaborated from time 
series of three-phase electric current measured from an electric arc furnace. Electric arc furnaces are 
industrial plants used for the production of steel by a melting process of metallic scraps. The V-I 
characteristic of the EAF under analysis is the result of a nonlinear process having memory. In fact, the 
value of the electrical current is determined as a function of both voltage and internal state, which in 
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turn is determined by considering the history of the states of the system. Such a process gives rise to an 
electrical phenomenon called hysteresis. 

We initially filtered the original signals via a Fourier based detrending, which we have performed in 
order to remove the implicit sinusoidal trends due to the alternate current. Our analysis showed that the 
memory effects of EAF process are suitably encoded into the filtered time series in the form of long-terms, 
positive correlations, which are quantified by the high values found for the Hurst exponent. Successively, 
we have shown the presence of a fair multifractal signature in all three time series, confirming thus the 
complexity of such noisy signals. Finally, shuffling of the (filtered) time series demonstrated that the 
multifractal signature was entirely due to the long-term correlations of the signals. 

Future research works will focus on the study of the long-term and multifractal properties of the 
cross-correlations between current and voltage time series related to the electric arc furnace. 

Multifractal detrended fluctuation analysis 

The MF-DFA procedure is extensively described in Ref. [19]. In this work, we partially relied on 
the software provided by Ihlen US). The method consists of five steps, three of which are identical to 
the DFA version. Given a time series X). of length N with compact support, the MF-DFA procedure is 
performed with the following steps: 

• Step 1 : The series Y{i) is computed as the cumulative sum (profile) of the series Xk'- 

i 

Yii) = '^[xk- {x}], i = (1) 

k=l 


• Step 2 : The series Y{i) is divided in Ng = mt{N/s) non-overlapping segments of equal length s. 
To account for the possible non-divisibility of N by s, the operation is repeated in reverse order by 
starting from the opposite end of the series, thus obtaining a total of 2Ns segments. 

• Step 3 : The local detrending operation is executed by performing a suitable polynomial fitting on 
each of the 2Ns segments. Then the variance is determined as 

s) = g ^ |^[(^ - l)s + *] - : (2) 

for each segment v = 1,..., Ng and 

s) = ^ |r[iV - (:/ - Ng)s + z] - 2/4 z)| (3) 

for V = Ng -|- 1,..., 2fVs, where y^{i) is the fitted polynomial in segment v. The order m of the 
fitting polynomial, y,y{i), determines the capability of the (MF-)DFA in eliminating trends in the 
series, thus it has to be tuned according to the expected maximum trending order of the time 
series. 


• Step 4 ■ The gth-order average of the variance over all segments is evaluated as 




2Ns 

E 


19/2 


1/9 


(4) 


with g G K. The g-dependence of the fluctuations function Fq{s) allows to highlight the contribu¬ 
tions of fluctuations at different magnitude orders. For q > 0 only the larger fluctuations contribute 
mostly to the average in Eq. conversely, for g < 0 the magnitude of the smaller fluctuations is 
enhanced. For q = 2 the standard DFA procedure is obtained. The case q = 0 cannot be computed 
with the averaging form in Eq. [^ and so a logarithmic form has to be employed. 


Fo{s) = exp 


1 


2N, 


^In 


I/=l 


(5) 


The steps 2 to 4 have to be repeated for different time scales s, where all values of s have to be 
chosen such that s > m + 2 to allow for a meaningful fitting of data. 


7 



• Step 5 : The scaling behaviour of the fluctuation functions can be determined by analyzing log-log 
plots of Fq{s) versus s for each value of q. If the series xi is long-range power-law correlated, Fq{s) 
is approximated (for large values of s) by the form 

( 6 ) 

The exponent H{q) is the generalized Hurst exponent; for q = 2 it reduces to the standard Hurst 
exponent, iJ. When the considered time series is monofractal, i.e., it shows a uniform scaling over 
all magnitude scales of the fluctuations, F[(q) is independent of q. On the contrary, when the small 
fluctuations scale differently from the large ones, the dependency of H{q) on q becomes apparent and 
the series can be considered multifractal. 

Starting from Eq. and using Eq. it is straightforward to obtain 

N/s 

i/=i 

where, for simplicity, it has been assumed that the length N of the series is a multiple of the scale s, 
such that Ns = N/s. The exponent 

r(g) = qHiq) - 1 (8) 

is the q-order mass exponent (also called Renyi scaling exponent) of the generalized partition function, 
Zq{s). It is worth mentioning that t(2) gives an estimation of the correlation dimension. The multifractal 
spectrum, denoted as /{■), provides a compact characterization of the multifractality of the time series. 
Such a function is obtained via the Legendre transform of 'r(q), 

f{a) =qa-T{q), (9) 

where a, called singularity exponent, is equal to the derivative r'(q). Using Eq. §it is possible to directly 
relate a and f{a) to F[{q), obtaining: 

a = F[{q) + qF['(q) and f{a)=q[a-F[{q)] + l. (10) 

The multifractal spectrum Q allows to infer important information regarding the degree of multi¬ 
fractality and the specific sensitivity of the time series to fluctuations of high/low magnitudes. In fact, 
the width of the support of /(•), defined as 

Aa = a{q-) - a{q+), (11) 

is a direct and important quantitative indicator of the multifractal signature in the series. 
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